A precision functional atlas of personalized network topography and probabilities.

Although the general location of functional neural networks is similar across individuals, there is vast person-to-person topographic variability. To capture this, we implemented precision brain mapping functional magnetic resonance imaging methods to establish an open-source, method-flexible set of precision functional network atlases-the Masonic Institute for the Developing Brain (MIDB) Precision Brain Atlas. This atlas is an evolving resource comprising 53,273 individual-specific network maps, from more than 9,900 individuals, across ages and cohorts, including the Adolescent Brain Cognitive Development study, the Developmental Human Connectome Project and others. We also generated probabilistic network maps across multiple ages and integration zones (using a new overlapping mapping technique, Overlapping MultiNetwork Imaging). Using regions of high network invariance improved the reproducibility of executive function statistical maps in brain-wide associations compared to group average-based parcellations. Finally, we provide a potential use case for probabilistic maps for targeted neuromodulation. The atlas is expandable to alternative datasets with an online interface encouraging the scientific community to explore and contribute to understanding the human brain function more precisely.

Multi-level fMRI analysis applied to hemispheric specialization in the language network, functional areas, and their behavioral correlations in the ABCD sample.

Prior research suggests that the organization of the language network in the brain is left-dominant and becomes more lateralized with age and increasing language skill. The age at which specific components of the language network become adult-like varies depending on the abilities they subserve. So far, a large, developmental study has not included a language task paradigm, so we introduce a method to study resting-state laterality in the Adolescent Brain Cognitive Development (ABCD) study. Our approach mixes source timeseries between left and right homotopes of the (1) inferior frontal and (2) middle temporal gyri and (3) a region we term "Wernicke's area" near the supramarginal gyrus. Our large subset sample size of ABCD (n = 6153) allows improved reliability and validity compared to previous, smaller studies of brain-behavior associations. We show that behavioral metrics from the NIH Youth Toolbox and other resources are differentially related to tasks with a larger linguistic component over ones with less (e.g., executive function-dominant tasks). These baseline characteristics of hemispheric specialization in youth are critical for future work determining the correspondence of lateralization with language onset in earlier stages of development.

Cumulative Effects of Resting-State Connectivity Across All Brain Networks Significantly Correlate with Attention-Deficit Hyperactivity Disorder Symptoms.

Identification of replicable neuroimaging correlates of attention-deficit hyperactivity disorder (ADHD) has been hindered by small sample sizes, small effects, and heterogeneity of methods. Given evidence that ADHD is associated with alterations in widely distributed brain networks and the small effects of individual brain features, a whole-brain perspective focusing on cumulative effects is warranted. The use of large, multisite samples is crucial for improving reproducibility and clinical utility of brain-wide MRI association studies. To address this, a polyneuro risk score (PNRS) representing cumulative, brain-wide, ADHD-associated resting-state functional connectivity was constructed and validated using data from the Adolescent Brain Cognitive Development (ABCD, N = 5,543, 51.5% female) study, and was further tested in the independent Oregon-ADHD-1000 case-control cohort (N = 553, 37.4% female). The ADHD PNRS was significantly associated with ADHD symptoms in both cohorts after accounting for relevant covariates (p < 0.001). The most predictive PNRS involved all brain networks, though the strongest effects were concentrated among the default mode and cingulo-opercular networks. In the longitudinal Oregon-ADHD-1000, non-ADHD youth had significantly lower PNRS (Cohen's d = -0.318, robust p = 5.5 × 10-4) than those with persistent ADHD (age 7-19). The PNRS, however, did not mediate polygenic risk for ADHD. Brain-wide connectivity was robustly associated with ADHD symptoms in two independent cohorts, providing further evidence of widespread dysconnectivity in ADHD. Evaluation in enriched samples demonstrates the promise of the PNRS approach for improving reproducibility in neuroimaging studies and unraveling the complex relationships between brain connectivity and behavioral disorders.

Heritability estimation of subcortical volumes in a multi-ethnic multi-site cohort study.

Heritability of regional subcortical brain volumes (rSBVs) describes the role of genetics in middle and inner brain development. rSBVs are highly heritable in adults but are not characterized well in adolescents. The Adolescent Brain Cognitive Development study (ABCD), taken over 22 US sites, provides data to characterize the heritability of subcortical structures in adolescence. In ABCD, site-specific effects co-occur with genetic effects which can bias heritability estimates. Existing methods adjusting for site effects require additional steps to adjust for site effects and can lead to inconsistent estimation. We propose a random-effect model-based method of moments approach that is a single step estimator and is a theoretically consistent estimator even when sites are imbalanced and performs well under simulations. We compare methods on rSBVs from ABCD. The proposed approach yielded heritability estimates similar to previous results derived from single-site studies. The cerebellum cortex and hippocampus were the most heritable regions (> 50%).

Neural Circuit Markers of Familial Risk for Depression Among Healthy Youths in the Adolescent Brain Cognitive Development Study.

BACKGROUND: Family history of depression is a robust predictor of early-onset depression, which may confer risk through alterations in neural circuits that have been implicated in reward and emotional processing. These alterations may be evident in youths who are at familial risk for depression but who do not currently have depression. However, the identification of robust and replicable findings has been hindered by few studies and small sample sizes. In the current study, we sought to identify functional connectivity (FC) patterns associated with familial risk for depression. METHODS: Participants included healthy (i.e., no lifetime psychiatric diagnoses) youths at high familial risk for depression (HR) (n = 754; at least one parent with a history of depression) and healthy youths at low familial risk for psychiatric problems (LR) (n = 1745; no parental history of psychopathology) who were 9 to 10 years of age and from the Adolescent Brain Cognitive Development (ABCD) Study sample. We conducted whole-brain seed-to-voxel analyses to examine group differences in resting-state FC with the amygdala, caudate, nucleus accumbens, and putamen. We hypothesized that HR youths would exhibit global amygdala hyperconnectivity and striatal hypoconnectivity patterns primarily driven by maternal risk. RESULTS: HR youths exhibited weaker caudate-angular gyrus FC than LR youths (α = 0.04, Cohen's d = 0.17). HR youths with a history of maternal depression specifically exhibited weaker caudate-angular gyrus FC (α = 0.03, Cohen's d = 0.19) as well as weaker caudate-dorsolateral prefrontal cortex FC (α = 0.04, Cohen's d = 0.21) than LR youths. CONCLUSIONS: Weaker striatal connectivity may be related to heightened familial risk for depression, primarily driven by maternal history. Identifying brain-based markers of depression risk in youths can inform approaches to improving early detection, diagnosis, and treatment.

Personalized functional brain network topography is associated with individual differences in youth cognition.

Individual differences in cognition during childhood are associated with important social, physical, and mental health outcomes in adolescence and adulthood. Given that cortical surface arealization during development reflects the brain's functional prioritization, quantifying variation in the topography of functional brain networks across the developing cortex may provide insight regarding individual differences in cognition. We test this idea by defining personalized functional networks (PFNs) that account for interindividual heterogeneity in functional brain network topography in 9-10 year olds from the Adolescent Brain Cognitive Development℠ Study. Across matched discovery (n = 3525) and replication (n = 3447) samples, the total cortical representation of fronto-parietal PFNs positively correlates with general cognition. Cross-validated ridge regressions trained on PFN topography predict cognition in unseen data across domains, with prediction accuracy increasing along the cortex's sensorimotor-association organizational axis. These results establish that functional network topography heterogeneity is associated with individual differences in cognition before the critical transition into adolescence.

Early-in-life isoflurane exposure alters resting-state functional connectivity in juvenile non-human primates.

BACKGROUND: Clinical studies suggest that anaesthesia exposure early in life affects neurobehavioural development. We designed a non-human primate (NHP) study to evaluate cognitive, behavioural, and brain functional and structural alterations after isoflurane exposure during infancy. These NHPs displayed decreased close social behaviour and increased astrogliosis in specific brain regions, most notably in the amygdala. Here we hypothesise that resting-state functional connectivity MRI can detect alterations in connectivity of brain areas that relate to these social behaviours and astrogliosis. METHODS: Imaging was performed in 2-yr-old NHPs under light anaesthesia, after early-in-life (postnatal days 6-12) exposure to 5 h of isoflurane either one or three times, or to room air. Brain images were segmented into 82 regions of interest; the amygdala and the posterior cingulate cortex were chosen for a seed-based resting-state functional connectivity MRI analysis. RESULTS: We found differences between groups in resting-state functional connectivity of the amygdala and the auditory cortices, medial premotor cortex, and posterior cingulate cortex. There were also alterations in resting-state functional connectivity between the posterior cingulate cortex and secondary auditory, polar prefrontal, and temporal cortices, and the anterior insula. Relationships were identified between resting-state functional connectivity alterations and the decrease in close social behaviour and increased astrogliosis. CONCLUSIONS: Early-in-life anaesthesia exposure in NHPs is associated with resting-state functional connectivity alterations of the amygdala and the posterior cingulate cortex with other brain regions, evident at the juvenile age of 2 yr. These changes in resting-state functional connectivity correlate with the decrease in close social behaviour and increased astrogliosis. Using resting-state functional connectivity MRI to study the neuronal underpinnings of early-in-life anaesthesia-induced behavioural alterations could facilitate development of a biomarker for anaesthesia-induced developmental neurotoxicity.

Using synthetic MR images for distortion correction.

Functional MRI (fMRI) data acquired using echo-planar imaging (EPI) are highly distorted by magnetic field inhomogeneities. Distortion and differences in image contrast between EPI and T1-weighted and T2-weighted (T1w/T2w) images makes their alignment a challenge. Typically, field map data are used to correct EPI distortions. Alignments achieved with field maps can vary greatly and depends on the quality of field map data. However, many public datasets lack field map data entirely. Additionally, reliable field map data is often difficult to acquire in high-motion pediatric or developmental cohorts. To address this, we developed Synth, a software package for distortion correction and cross-modal image registration that does not require field map data. Synth combines information from T1w and T2w anatomical images to construct an idealized undistorted synthetic image with similar contrast properties to EPI data. This synthetic image acts as an effective reference for individual-specific distortion correction. Using pediatric (ABCD: Adolescent Brain Cognitive Development) and adult (MSC: Midnight Scan Club; HCP: Human Connectome Project) data, we demonstrate that Synth performs comparably to field map distortion correction approaches, and often outperforms them. Field map-less distortion correction with Synth allows accurate and precise registration of fMRI data with missing or corrupted field map information.

BIBSNet: A Deep Learning Baby Image Brain Segmentation Network for MRI Scans.

Objectives: Brain segmentation of infant magnetic resonance (MR) images is vitally important in studying developmental mental health and disease. The infant brain undergoes many changes throughout the first years of postnatal life, making tissue segmentation difficult for most existing algorithms. Here, we introduce a deep neural network BIBSNet (Baby and Infant Brain Segmentation Neural Network), an open-source, community-driven model that relies on data augmentation and a large sample size of manually annotated images to facilitate the production of robust and generalizable brain segmentations. Experimental Design: Included in model training and testing were MR brain images on 84 participants with an age range of 0-8 months (median postmenstrual ages of 13.57 months). Using manually annotated real and synthetic segmentation images, the model was trained using a 10-fold cross-validation procedure. Testing occurred on MRI data processed with the DCAN labs infant-ABCD-BIDS processing pipeline using segmentations produced from gold standard manual annotation, joint-label fusion (JLF), and BIBSNet to assess model performance. Principal Observations: Using group analyses, results suggest that cortical metrics produced using BIBSNet segmentations outperforms JLF segmentations. Additionally, when analyzing individual differences, BIBSNet segmentations perform even better. Conclusions: BIBSNet segmentation shows marked improvement over JLF segmentations across all age groups analyzed. The BIBSNet model is 600x faster compared to JLF and can be easily included in other processing pipelines.

Polyneuro risk scores capture widely distributed connectivity patterns of cognition.

Resting-state functional connectivity (RSFC) is a powerful tool for characterizing brain changes, but it has yet to reliably predict higher-order cognition. This may be attributed to small effect sizes of such brain-behavior relationships, which can lead to underpowered, variable results when utilizing typical sample sizes (N∼25). Inspired by techniques in genomics, we implement the polyneuro risk score (PNRS) framework - the application of multivariate techniques to RSFC data and validation in an independent sample. Utilizing the Adolescent Brain Cognitive Development® cohort split into two datasets, we explore the framework's ability to reliably capture brain-behavior relationships across 3 cognitive scores - general ability, executive function, learning & memory. The weight and significance of each connection is assessed in the first dataset, and a PNRS is calculated for each participant in the second. Results support the PNRS framework as a suitable methodology to inspect the distribution of connections contributing towards behavior, with explained variance ranging from 1.0 % to 21.4 %. For the outcomes assessed, the framework reveals globally distributed, rather than localized, patterns of predictive connections. Larger samples are likely necessary to systematically identify the specific connections contributing towards complex outcomes. The PNRS framework could be applied translationally to identify neurologically distinct subtypes of neurodevelopmental disorders.

Carotenoids improve the development of cerebral cortical networks in formula-fed infant macaques.

Nutrition during the first years of life has a significant impact on brain development. This study characterized differences in brain maturation from birth to 6 months of life in infant macaques fed formulas differing in content of lutein, ß-carotene, and other carotenoids using Magnetic Resonance Imaging to measure functional connectivity. We observed differences in functional connectivity based on the interaction of diet, age and brain networks. Post hoc analysis revealed significant diet-specific differences between insular-opercular and somatomotor networks at 2 months of age, dorsal attention and somatomotor at 4 months of age, and within somatomotor and between somatomotor-visual and auditory-dorsal attention networks at 6 months of age. Overall, we found a larger divergence in connectivity from the breastfeeding group in infant macaques fed formula containing no supplemental carotenoids in comparison to those fed formula supplemented with carotenoids. These findings suggest that carotenoid formula supplementation influences functional brain development.

Curation of BIDS (CuBIDS): A workflow and software package for streamlining reproducible curation of large BIDS datasets.

The Brain Imaging Data Structure (BIDS) is a specification accompanied by a software ecosystem that was designed to create reproducible and automated workflows for processing neuroimaging data. BIDS Apps flexibly build workflows based on the metadata detected in a dataset. However, even BIDS valid metadata can include incorrect values or omissions that result in inconsistent processing across sessions. Additionally, in large-scale, heterogeneous neuroimaging datasets, hidden variability in metadata is difficult to detect and classify. To address these challenges, we created a Python-based software package titled "Curation of BIDS" (CuBIDS), which provides an intuitive workflow that helps users validate and manage the curation of their neuroimaging datasets. CuBIDS includes a robust implementation of BIDS validation that scales to large samples and incorporates DataLad--a version control software package for data--as an optional dependency to ensure reproducibility and provenance tracking throughout the entire curation process. CuBIDS provides tools to help users perform quality control on their images' metadata and identify unique combinations of imaging parameters. Users can then execute BIDS Apps on a subset of participants that represent the full range of acquisition parameters that are present, accelerating pipeline testing on large datasets.

An open-access accelerated adult equivalent of the ABCD Study neuroimaging dataset (a-ABCD).

As public access to longitudinal developmental datasets like the Adolescent Brain Cognitive Development StudySM (ABCD Study®) increases, so too does the need for resources to benchmark time-dependent effects. Scan-to-scan changes observed with repeated imaging may reflect development but may also reflect practice effects, day-to-day variability in psychological states, and/or measurement noise. Resources that allow disentangling these time-dependent effects will be useful in quantifying actual developmental change. We present an accelerated adult equivalent of the ABCD Study dataset (a-ABCD) using an identical imaging protocol to acquire magnetic resonance imaging (MRI) structural, diffusion-weighted, resting-state and task-based data from eight adults scanned five times over five weeks. We report on the task-based imaging data (n = 7). In-scanner stop-signal (SST), monetary incentive delay (MID), and emotional n-back (EN-back) task behavioral performance did not change across sessions. Post-scan recognition memory for emotional n-back stimuli, however, did improve as participants became more familiar with the stimuli. Functional MRI analyses revealed that patterns of task-based activation reflecting inhibitory control in the SST, reward success in the MID task, and working memory in the EN-back task were more similar within individuals across repeated scan sessions than between individuals. Within-subject, activity was more consistent across sessions during the EN-back task than in the SST and MID task, demonstrating differences in fMRI data reliability as a function of task. The a-ABCD dataset provides a unique testbed for characterizing the reliability of brain function, structure, and behavior across imaging modalities in adulthood and benchmarking neurodevelopmental change observed in the open-access ABCD Study.

Synthesizing pseudo-T2w images to recapture missing data in neonatal neuroimaging with applications in rs-fMRI.

T1- and T2-weighted (T1w and T2w) images are essential for tissue classification and anatomical localization in Magnetic Resonance Imaging (MRI) analyses. However, these anatomical data can be challenging to acquire in non-sedated neonatal cohorts, which are prone to high amplitude movement and display lower tissue contrast than adults. As a result, one of these modalities may be missing or of such poor quality that they cannot be used for accurate image processing, resulting in subject loss. While recent literature attempts to overcome these issues in adult populations using synthetic imaging approaches, evaluation of the efficacy of these methods in pediatric populations and the impact of these techniques in conventional MR analyses has not been performed. In this work, we present two novel methods to generate pseudo-T2w images: the first is based in deep learning and expands upon previous models to 3D imaging without the requirement of paired data, the second is based in nonlinear multi-atlas registration providing a computationally lightweight alternative. We demonstrate the anatomical accuracy of pseudo-T2w images and their efficacy in existing MR processing pipelines in two independent neonatal cohorts. Critically, we show that implementing these pseudo-T2w methods in resting-state functional MRI analyses produces virtually identical functional connectivity results when compared to those resulting from T2w images, confirming their utility in infant MRI studies for salvaging otherwise lost subject data.

Reproducible brain-wide association studies require thousands of individuals.

Magnetic resonance imaging (MRI) has transformed our understanding of the human brain through well-replicated mapping of abilities to specific structures (for example, lesion studies) and functions1-3 (for example, task functional MRI (fMRI)). Mental health research and care have yet to realize similar advances from MRI. A primary challenge has been replicating associations between inter-individual differences in brain structure or function and complex cognitive or mental health phenotypes (brain-wide association studies (BWAS)). Such BWAS have typically relied on sample sizes appropriate for classical brain mapping4 (the median neuroimaging study sample size is about 25), but potentially too small for capturing reproducible brain-behavioural phenotype associations5,6. Here we used three of the largest neuroimaging datasets currently available-with a total sample size of around 50,000 individuals-to quantify BWAS effect sizes and reproducibility as a function of sample size. BWAS associations were smaller than previously thought, resulting in statistically underpowered studies, inflated effect sizes and replication failures at typical sample sizes. As sample sizes grew into the thousands, replication rates began to improve and effect size inflation decreased. More robust BWAS effects were detected for functional MRI (versus structural), cognitive tests (versus mental health questionnaires) and multivariate methods (versus univariate). Smaller than expected brain-phenotype associations and variability across population subsamples can explain widespread BWAS replication failures. In contrast to non-BWAS approaches with larger effects (for example, lesions, interventions and within-person), BWAS reproducibility requires samples with thousands of individuals.

Filtering respiratory motion artifact from resting state fMRI data in infant and toddler populations.

The importance of motion correction when processing resting state functional magnetic resonance imaging (rs-fMRI) data is well-established in adult cohorts. This includes adjustments based on self-limited, large amplitude subject head motion, as well as factitious rhythmic motion induced by respiration. In adults, such respiration artifact can be effectively removed by applying a notch filter to the motion trace, resulting in higher amounts of data retained after frame censoring (e.g., "scrubbing") and more reliable correlation values. Due to the unique physiological and behavioral characteristics of infants and toddlers, rs-fMRI processing pipelines, including methods to identify and remove colored noise due to subject motion, must be appropriately modified to accurately reflect true neuronal signal. These younger cohorts are characterized by higher respiration rates and lower-amplitude head movements than adults; thus, the presence and significance of comparable respiratory artifact and the subsequent necessity of applying similar techniques remain unknown. Herein, we identify and characterize the consistent presence of respiratory artifact in rs-fMRI data collected during natural sleep in infants and toddlers across two independent cohorts (aged 8-24 months) analyzed using different pipelines. We further demonstrate how removing this artifact using an age-specific notch filter allows for both improved data quality and data retention in measured results. Importantly, this work reveals the critical need to identify and address respiratory-driven head motion in fMRI data acquired in young populations through the use of age-specific motion filters as a mechanism to optimize the accuracy of measured results in this population.

QSIPrep: an integrative platform for preprocessing and reconstructing diffusion MRI data.

Diffusion-weighted magnetic resonance imaging (dMRI) is the primary method for noninvasively studying the organization of white matter in the human brain. Here we introduce QSIPrep, an integrative software platform for the processing of diffusion images that is compatible with nearly all dMRI sampling schemes. Drawing on a diverse set of software suites to capitalize on their complementary strengths, QSIPrep facilitates the implementation of best practices for processing of diffusion images.

Characterizing the impact of adversity, abuse, and neglect on adolescent amygdala resting-state functional connectivity.

Characterizing typologies of childhood adversity may inform the development of risk profiles and corresponding interventions aimed at mitigating its lifelong consequences. A neurobiological grounding of these typologies requires systematic comparisons of neural structure and function among individuals with different exposure histories. Using seed-to-whole brain analyses, this study examined associations between childhood adversity and amygdala resting-state functional connectivity (rs-fc) in adolescents aged 11-19 years across three independent studies (N = 223; 127 adversity group) in both general and dimensional models of adversity (comparing abuse and neglect). In a general model, adversity was associated with altered amygdala rs-fc with clusters within the left anterior lateral prefrontal cortex. In a dimensional model, abuse was associated with altered amygdala rs-fc within the orbitofrontal cortex, dorsal precuneus, posterior cingulate cortex, and dorsal anterior cingulate cortex/anterior mid-cingulate cortex, as well as within the dorsal attention, visual, and somatomotor networks. Neglect was associated with altered amygdala rs-fc with the hippocampus, supplementary motor cortex, temporoparietal junction, and regions within the dorsal attention network. Both general and dimensional models revealed unique regions, potentially reflecting pathways by which distinct histories of adversity may influence adolescent behavior, cognition, and psychopathology.

Heterogeneity of executive function revealed by a functional random forest approach across ADHD and ASD.

BACKGROUND: Those with autism spectrum disorder (ASD) and/or attention-deficit-hyperactivity disorder (ADHD) exhibit symptoms of hyperactivity and inattention, causing significant hardships for families and society. A potential mechanism involved in these conditions is atypical executive function (EF). Inconsistent findings highlight that EF features may be shared or distinct across ADHD and ASD. With ADHD and ASD each also being heterogeneous, we hypothesized that there may be nested subgroups across disorders with shared or unique underlying mechanisms. METHODS: Participants (N = 130) included adolescents aged 7-16 with ASD (n = 64) and ADHD (n = 66). Typically developing (TD) participants (n = 28) were included for a comparative secondary sub-group analysis. Parents completed the K-SADS and youth completed an extended battery of executive and other cognitive measures. A two stage hybrid machine learning tool called functional random forest (FRF) was applied as a classification approach and then subsequently to subgroup identification. We input 43 EF variables to the classification step, a supervised random forest procedure in which the features estimated either hyperactive or inattentive ADHD symptoms per model. The FRF then produced proximity matrices and identified optimal subgroups via the infomap algorithm (a type of community detection derived from graph theory). Resting state functional connectivity MRI (rs-fMRI) was used to evaluate the neurobiological validity of the resulting subgroups. RESULTS: Both hyperactive (Mean absolute error (MAE) = 0.72, Null model MAE = 0.8826, (t(58) = -4.9, p < .001) and inattentive (MAE = 0.7, Null model MAE = 0.85, t(58) = -4.4, p < .001) symptoms were predicted better than chance by the EF features selected. Subgroup identification was robust (Hyperactive: Q = 0.2356, p < .001; Inattentive: Q = 0.2350, p < .001). Two subgroups representing severe and mild symptomology were identified for each symptom domain. Neuroimaging data revealed that the subgroups and TD participants significantly differed within and between multiple functional brain networks, but no consistent "severity" patterns of over or under connectivity were observed between subgroups and TD. CONCLUSION: The FRF estimated hyperactive/inattentive symptoms and identified 2 distinct subgroups per model, revealing distinct neurocognitive profiles of Severe and Mild EF performance per model. Differences in functional connectivity between subgroups did not appear to follow a severity pattern based on symptom expression, suggesting a more complex mechanistic interaction that cannot be attributed to symptom presentation alone.